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BACKGROUND OF THE INVENTION 
1. Field of the Invention 

The invention relates generally to the field of orthodontics and, more 
10 particularly, to computer-automated development of an orthodontic treatment plan and 
appliance. 

Repositioning teeth for aesthetic or other reasons is accomplished 
conventionally by wearing what are commonly referred to as "braces." Braces comprise a 
variety of appliances such as brackets, archwires, ligatures, and O-rings. Attaching the 
1 5 appliances to a patient's teeth is a tedious and time consuming enterprise requiring many 
meetings with the treating orthodontist. Consequently, conventional orthodontic 
treatment limits an orthodontist's patient capacity and makes orthodontic treatment quite 
expensive. 

Before fastening braces to a patient's teeth, at least one appointment is 
20 typically scheduled with the orthodontist, dentist, and/or X-ray laboratory so that X-rays 
and photographs of the patient's teeth and jaw structure can be taken. Also during this 
preliminary meeting, or possibly at a later meeting, an alginate mold of the patient's teeth 
is typically made. This mold provides a model of the patient's teeth that the orthodontist 
uses in conjunction with the X-rays and photographs to formulate a treatment strategy. 
25 The orthodontist then typically schedules one or more appointments during which braces 
will be attached to the patient's teeth. 

At the meeting during which braces are first attached, the teeth surfaces 
are initially treated with a weak acid. The acid optimizes the adhesion properties of the 
teeth surfaces for brackets and bands that are to be bonded to them. The brackets and 
30 bands serve as anchors for other appliances to be added later. After the acid step, the 

brackets and bands are cemented to the patient's teeth using a suitable bonding material. 
No force-inducing appliances are added until the cement is set. For this reason, it is 




common for the orthodontist to schedule a later appointment to ensure that the brackets 
and bands are well bonded to the teeth. 

The primary force-inducing appliance in a conventional set of braces is the 
archwire. The archwire is resilient and is attached to the brackets by way of slots in the 
5 brackets. The archwire links the brackets together and exerts forces on them to move the 
teeth over time. Twisted wires or elastomeric O-rings are commonly used to reinforce 
attachment of the archwire to the brackets. Attachment of the archwire to the brackets is 
known in the art of orthodontia as "ligation" and wires used in this procedure are called 
"ligatures." The elastomeric O-rings are called "plastics." 

1 0 After the archwire is in place, periodic meetings with the orthodontist are 

required, during which the patient's braces will be adjusted by installing a different 
archwire having different force-inducing properties or by replacing or tightening existing 
ligatures. Typically, these meetings are scheduled every three to six weeks. 

As the above illustrates, the use of conventional braces is a tedious and 

15 time consuming process and requires many visits to the orthodontist's office. Moreover, 
from the patient's perspective, the use of braces is unsightly, uncomfortable, presents a 
risk of infection, and makes brushing, flossing, and other dental hygiene procedures 
difficult. 

For these reasons, it would be desirable to provide alternative methods and 
20 systems for repositioning teeth. Such methods and systems should be economical, and in 
particular should reduce the amount of time required by the orthodontist in planning and 
overseeing each individual patient. The methods and systems should also be more 
acceptable to the patient, in particular being less visible, less uncomfortable, less prone to 
infection, and more compatible with daily dental hygiene. At least some of these 
25 objectives will be met by the methods and systems of the present invention described 
hereinafter. 

2. Description of the Background Art 

Tooth positioners for finishing orthodontic treatment are described by 
Kesling in the Am. J. Orthod. Oral. Surg. 57:297-304 (1945) and 52:285-293 (1946). The 
30 use of silicone positioners for the comprehensive orthodontic realignment of a patient's 
teeth is described in Warunek et al. (1989) J. Clin. Orthod. 25:694-700. Clear plastic 
retainers for finishing and maintaining tooth positions are commercially available from 
Raintree Essix, Inc., New Orleans, Louisiana 70125, and Tru-Tain Plastics, Rochester, 
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Minnesota 55902. The manufacture of orthodontic positioners is described in U.S. Patent 
Nos. 5,186,623; 5,059,1 18; 5,055,039; 5,035,613; 4,856,991; 4,798,534; and 4,755,139. 

Other publications describing the fabrication and use of dental positioners 
include Kleemann and Janssen (1996) J. Clin. Orthodon. 30:673-680; Cureton (1996) J. 
5 Clin. Orthodon. 30:390-395; Chiappone (1980)/. Clin. Orthodon. 7*121-133; Shilliday 
(1971) Am. J. Orthodontics 59:596-599; Wells (1970) Am. J. Orthodontics 55:351-366; 
and Cottingham (1969) Am. J. Orthodontics 55:23-31. 

Kuroda et al. (1996) Am. J. Orthodontics 770:365-369 describes a method 
for laser scanning a plaster dental cast to produce a digital image of the cast. See also 
10 U.S. Patent No. 5,605,459. 

U.S. Patent Nos. 5,533,895; 5,474,448; 5,454,717; 5,447,432; 5,431,562; 
5,395,238; 5,368,478; and 5,139.419, assigned to Ormco Corporation, describe methods 
for manipulating digital images of teeth for designing orthodontic appliances. 

U.S. Patent No. 5,01 1,405 describes a method for digitally imaging a tooth 
1 5 and determining optimum bracket positioning for orthodontic treatment. Laser scanning 
of a molded tooth to produce a three-dimensional model is described in U.S. Patent No. 
5,338,198. U.S. Patent NO. 5,452,219 describes a method for laser scanning a tooth 
model and milling a tooth mold. Digital computer manipulation of tooth contours is 
described in U.S. Patent Nos. 5,607,305 and 5,587,912. Computerized digital imaging of 
20 the jaw is described in U.S. Patent Nos. 5,342,202 and 5,340,309. Other patents of 

interest include U.S. Patent Nos. 5,549,476; 5,382,164; 5,273,429; 4,936,862; 3,860,803; 
3,660,900; 5,645,421; 5,055,039; 4,798,534; 4,856,991; 5,035,613; 5,059,118; 5,186,623; 
and 4,755,139. 

SUMMARY OF THE INVENTION 

25 In one aspect, the invention relates to the computer-automated creation of 

a plan for repositioning an orthodontic patient's teeth. A computer receives an initial 
digital data set representing the patient's teeth at their initial positions and a final digital 
data set representing the teeth at their final positions. The computer uses the data sets to 
generate treatment paths along which the teeth will move from the initial positions to the 

30 final positions. 

In some implementations, the initial digital data set includes data obtained 
by scanning a physical model of the patient's teeth, such as by scanning a positive 
impression or a negative impression of the patient's teeth with a laser scanner or a 
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destructive scanner. The positive and negative impression may be scanned while 
interlocked with each other to provide more accurate data. The initial digital data set also 
may include volume image data of the patient's teeth, which the computer can convert 
into a 3D geometric model of the tooth surfaces, for example using a conventional 
5 marching cubes technique. The computer also can be used to segment the initial digital 
data set automatically into individual tooth models, such as by performing a feature 
detection or matching operation on the image data. In some embodiments, the individual 
tooth models include data representing hidden tooth surfaces, such as roots imaged 
through x-ray, CT scan, or MRI techniques. Tooth roots and hidden surfaces also can be 

10 extrapolated from the visible surfaces of the patient's teeth. 

In other embodiments, the computer applies a set of rules to detect 
collisions that will occur as the patient's teeth move along the treatment paths. One 
technique for collision detection is the creation of a collision buffer between two teeth at 
a given step along the treatment path. The computer also can be used to detect improper 

1 5 bite occlusions that will occur as the patient's teeth move along the treatment paths. 
Other embodiments allow the computer to render a three-dimensional (3D) graphical 
representation of the teeth at any selected treatment step. The computer also can be used 
to animate the graphical representation of the teeth to provide a visual display of the 
movement of the teeth along the treatment paths. A VCR metaphor in the graphical user 

20 interface allows the user to control the animation. Level-of-detail compression can be 

used to improve the speed at with the 3D image of the teeth is rendered. Moreover, some 
embodiments allow the user to modify the underlying digital data set by repositioning a 
tooth in the 3D graphical representation. 

In another aspect, the invention involves generating three-dimensional 

25 models of individual teeth from an initial data set that contains a 3D representation of a 
group of teeth. A computer performs this task by identifying points in the initial data set 
corresponding to each individual tooth and then segmenting the initial data set into 
multiple data sets, each containing the points identified for one of the teeth. In some 
embodiments, the computer stores each data set as a 3D geometric model representing the 

30 visible surfaces of the corresponding tooth. The computer can be used to modify each 3D 
model to include hidden surfaces of the corresponding tooth. In other embodiments, the 
initial data set contains digital volume image data, and the computer converts the volume 
image data into a 3D geometric model by detecting volume elements in the image data 
between which a sharp transition in digital image value occurs. 



In another aspect, the invention relates to determining whether a patient's 



teeth can be moved from a first set of positions to a second set of positions. A computer 
performs this task by receiving a digital data set representing the teeth at the second set of 
positions and determining whether any of the teeth will collide while moving to the 
5 second set of positions. In some embodiments, the computer calculates distances between 
two of the teeth (a first tooth and a second tooth) by establishing a neutral projection 
plane between the first tooth and the second tooth, establishing a z-axis that is normal to 
the plane and that has a positive direction and a negative direction from each of a set of 
base points on the projection plane, computing a pair of signed distances comprising a 
1 0 first signed distance to the first tooth and a second signed distance to the second tooth, the 
signed distances being measured on a line passing through the base points and parallel to 
the z-axis, and determining that a collision will occur if any of the pair of signed distances 
indicates a collision. 



15 orthodontic patient's teeth. A computer receives a digital data set representing the teeth 
at recommended final positions, renders a three-dimensional (3D) graphical 
representation of the teeth at the recommended final positions, receives an instruction to 
reposition one of the teeth in response to a user's manipulation of the tooth in the 
graphical representation, and, in response to the instruction, modifies the digital data set 

20 to represent the teeth at the user-selected final positions. 



treatment plan for an orthodontic patient's teeth. A computer receives a digital data set 
representing the patient's upper teeth after treatment, receives a digital data set 
representing the patient's lower teeth after treatment, orients the data in the data sets to 
25 simulate the patient's bite occlusion, manipulates the data sets in a manner that simulates 
motion of human jaws, and detects collisions between the patient's upper teeth and lower 
teeth during the simulation of motion. The simulation of motion can be based on 
observed motion of typical human jaws or the patient's jaws. 



In another aspect, the invention relates to determining final positions for an 



In another aspect, the invention relates to analyzing a recommended 



BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1A illustrates a patient's jaw and provides a general indication of how 



teeth may be moved. 



FIG. IB illustrates a single tooth from FIG. 1 A and defines how tooth 



movement distances are determined. 
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FIG. 1C illustrates the jaw of FIG. 1A together with an incremental 
position adjustment appliance. 

FIG. 2 is a block diagram illustrating steps for producing a system of 
incremental position adjustment appliances. 
5 FIG. 3 is a block diagram setting forth the steps for manipulating an initial 

digital data set representing an initial tooth arrangement to produce a final digital data set 
corresponding to a desired final tooth arrangement. 

FIG. 4A is a flow chart illustrating an eraser tool for the methods herein. 

FIG. 4B illustrates the volume of space which is being erased by the 
10 program of FIG. 4A. 

FIG. 5 is a flow chart illustrating a program for matching high-resolution 
and low-resolution components in the manipulation of data sets of FIG. 3. 

FIG. 6A is a flow chart illustrating a program for performing the 
"detection" stage of the cusp detection algorithm. 
1 5 FIG. 6B is a flow chart illustrating a program for performing the 

"rejection" stage of the cusp detection algorithm. 

FIG. 7 illustrates a method for generating multiple intermediate digital 
data sets which are used for producing the adjustment appliances. 

FIG. 8A is a flow chart illustrating the steps performed by the path 
20 scheduling algorithm. 

FIG. 8B is a flow chart illustrating the steps for performing a "visibility" 

function. 

FIG. 8C is a flow chart illustrating the steps for performing a "children" 

function. 

25 FIG. 8D is a flow chart illustrating the steps for performing path 

scheduling step 128 of FIG. 8 A. 

FIG. 9A is a flow chart illustrating the steps for performing recursive 
collision testing during collision detection. 

FIG. 9B is a flow chart illustrating node splitting performed during 
30 collision detection. 

FIG. 9C is a flow chart illustrating steps for providing additional motion 
information to the collision detection process. 
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FIG. 1 0 illustrates alternative processes for producing a plurality of 
appliances utilizing digital data sets representing the intermediate and final appliance 
designs. 

FIG. 1 1 is a simplified block diagram of a data processing system. 
5 FIG. 12 is cross-sectional image of a plaster tooth casting in an epoxy 

mold. 

FIG. 13 is a flow chart of a process for forming one 3D image data set of 
teeth from two sets of image data. 

FIG. 14 is a flow chart of a process for forming a 3D surface mesh from 

10 3D image data. 

FIGS. 15 A, 15B 5 and 15C illustrate the positioning of teeth at various steps 
of an orthodontic treatment plan. 

FIG. 16 is a flow chart of a process for determining a tooth's path among 
intermediate positions during an orthodontic treatment plan. 
1 5 FIG. 1 7 is a flow chart of a process for optimizing the path of a tooth from 

an initial position to a final position during an orthodontic treatment plan. 

FIG. 18 is a diagram illustrating a buffering technique for use in a collision 
detection algorithm. 

FIG. 19 is a flow chart for a collision detection technique. 
20 FIG. 20 is a screen shot of a GUI display used to render 3D images of an 

orthodontic patient's teeth. 

FIGS. 21 A and 2 IB illustrate a technique for improving the downloading 
and rendering speed of an orthodontic image data file. 

DESCRIPTION OF THE SPECIFIC EMBODIMENTS 
25 Systems and methods are provided for moving teeth incrementally using a 

plurality of discrete appliances, where each appliance successively moves one or more of 
the patient's teeth by relatively small amounts. The tooth movements will be those 
normally associated with orthodontic treatment, including translation in all three 
orthogonal directions relative to a vertical centerline, rotation of the tooth centerline in the 
30 two orthodontic directions ("root angulation" and "torque"), as well as rotation about the 
centerline. 

Referring now to FIG. 1 A, a representative jaw 100 includes sixteen teeth, 
at least some of which are to be moved from an initial tooth arrangement to a final tooth 
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arrangement. To understand how the teeth may be moved, an arbitrary centerline (CL) is 
drawn through one of the teeth 102. With reference to this centerline (CL), the teeth may 
be moved in the orthogonal directions represented by axes 104, 106, and 108 (where 104 
is the centerline). The centerline may be rotated about the axis 1 08 (root angulation) and 
5 104 (torque) as indicated by arrows 1 10 and 1 12, respectively. Additionally, the tooth 
may be rotated about the centerline, as represented by arrow 1 1 4. Thus, all possible free- 
form motions of the tooth can be performed. 

Referring now to FIG. IB, the magnitude of any tooth movement is 
defined in terms of the maximum linear translation of any point P on a tooth 102. Each 

10 point Pj will undergo a cumulative translation as that tooth is moved in any of the 

orthogonal or rotational directions defined in FIG. 1 A. That is, while the point will 
usually follow a non-linear path, there will be a linear distance between any point in the 

tooth when determined at any two times during the treatment. Thus, an arbitrary point P \ 

may in fact undergo a true side- to-side translation as indicated by arrow di, while a 
1 5 second arbitrary point P 2 may travel along an arcuate path, resulting in a final translation 

d2- In many situations, the maximum permissible movement of a point Pj in any 

particular tooth is defined as the maximum linear translation of that point Pi on the tooth 
which undergoes the maximum movement for that tooth in any treatment step. 

One tool for a incrementally repositioning the teeth in a patient's jaw is a 

20 set of one or more adjustment appliances. Suitable appliances include any of the known 
positioners, retainers, or other removable appliances which are used for finishing and 
maintaining teeth positions in connection with conventional orthodontic treatment. As 
described below, a plurality of such appliances can be worn by a patient successively to 
achieve gradual tooth repositioning. A particularly advantageous appliance is the 

25 appliance 100, shown in FIG. 1C, which typically comprises a polymeric shell having a 
cavity shaped to receive and resiliently reposition teeth from one tooth arrangement to 
another tooth arrangement. The polymeric shell typically fits over all teeth present in the 
upper or lower jaw. Often, only some of the teeth will be repositioned while others will 
provide a base or anchor region for holding the repositioning appliance in place as it 

30 applies the resilient repositioning force against the tooth or teeth to be repositioned. In 
complex cases, however, many or most of the teeth will be repositioned at some point 
during the treatment. In such cases, the teeth which are moved can also serve as a base or 
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anchor region for holding the repositioning appliance. The gums and the palette also 
serve as an anchor region in some cases, thus allowing all or nearly all of the teeth to be 
repositioned simultaneously. 

The polymeric appliance 100 of FIG. 1C is preferably formed from a thin 
5 sheet of a suitable elastomeric polymeric, such as Tru-Tain 0.03 in. thermal forming 
dental material, marketed by Tru-Tain Plastics, Rochester, Minnesota 55902. In many 
cases, no wires or other means are provided for holding the appliance in place over the 
teeth. In some cases, however, it is necessary to provide individual attachments on the 
teeth with corresponding receptacles or apertures in the appliance 100 so that the 
10 appliance can apply forces that would not be possible or would be difficult to apply in 
the absence of such attachments. 

BUILDING A DIGITAL MODEL OF THE PATIENT'S TEETH 
Gathering data about the teeth 

Referring now to FIG. 2, a method for producing the incremental position 

1 5 adjustment appliances for subsequent use by a patient to reposition the patient's teeth will 
be described. As a first step, a digital data set representing an initial tooth arrangement is 
obtained, referred to hereinafter as the initial digital data set, or IDDS. The IDDS may be 
obtained in a variety of ways. For example, the patient's teeth may be scanned or imaged 
using well known technology, such as X-rays, three-dimensional X-rays, computer-aided 

20 tomographic images or data sets, and magnetic resonance images. Methods for digitizing 
such conventional images to produce useful data sets are well known and described in the 
patent and medical literature. Usually, however, a plaster cast of the patient's teeth is 
obtained by well known techniques, such as those described in Graber, Orthodontics: 
Principle and Practice, Second Edition, Saunders, Philadelphia, 1969, pp. 401-415. 

25 After the tooth casting is obtained, the casting is digitally scanned by a 

scanner, such as a non-contact type laser or destructive scanner or a contact-type scanner, 
to produce the IDDS. The data set produced by the scanner may be presented in any of a 
variety of digital formats to ensure compatibility with the software used to manipulate 
images represented by the data, as described in more detail below. General techniques for 

30 producing plaster casts of teeth and generating digital models using laser scanning 

techniques are described, for example, in U.S. Patent No. 5,605,459, the full disclosure of 
which is incorporated in this application by reference. 



9 




Suitable scanners include a variety of range acquisition systems, generally 
categorized by whether the acquisition process requires contact with the three 
dimensional object being scanned. Some contact-type scanners use probes having 
multiple degrees of translational and/or rotational freedom. A computer-readable (i.e., 
5 digital) representation of the sample object is generated by recording the physical 
displacement of the probe as it is drawn across the sample surface. 

Conventional non-contact-type scanners include reflective-type and 
transmissive-type systems. A wide variety of reflective systems are in use today, some of 
which utilize non-optical incident energy sources such as microwave radar or sonar. 
10 Others utilize optical energy. Non-contact-type systems that use reflected optical energy 
usually include special instrumentation that carry out certain measuring techniques (e.g., 
imaging radar, triangulation and interferometry). 

One type of non-contact scanner is an optical, reflective scanner, such as a 
laser scanner. Non-contact scanners such as this are inherently nondestructive (i.e., do 
15 not damage the sample object), generally are characterized by a relatively high capture 
resolution, and are capable of scanning a sample in a relatively short period of time. One 
such scanner is the Cyberware Model 15 scanner manufactured by Cyberware, Inc., 
Monterey, California. 

Both non-contact- type and contact- type scanners also can include color 
20 cameras which, when synchronized with the scanning capabilities, provide means for 

capturing, in digital format, color representations of the sample objects. The importance 
of this ability to capture not just the shape of the sample object but also its color is 
discussed below. 

Other scanners, such as the CSS- 1000 model destructive scanner produced 
25 by Capture Geometry Inside (CGI), Minneapolis, Minnesota, can provide more detailed 
and precise information about a patient's teeth than a typical range acquisition scanner 
can provide. In particular, a destructive scanner can image areas that are hidden or 
shielded from a range acquisition scanner and therefore may not be subject to adequate 
imaging. The CSS- 1000 scanner gathers image data for an object by repeatedly milling 
30 thin slices from the object and optically scanning the sequence of milled surfaces to create 
a sequence of 2D image slices, so none of the object's surfaces are hidden from the 
scanner. Image processing software combines the data from individual slices to form a 
data set representing the object, which later is converted into a digital representation of 
the surfaces of the object, as described below. 
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The destructive scanner may be used in conjunction with a laser scanner to 
create a digital model of a patient's teeth. For example, a laser scanner may be used first 
to build a low-resolution image of a patient's upper and lower arches coupled with the 
patient's wax bite, as described below. The destructive scanner then may be used to form 
5 high-resolution images of the individual arches. The data obtained by the laser scanner 
indicates the relation between the patient's upper and lower teeth, which later can be used 
to relate to each other the images generated by the destructive scanner and the digital 
models derived from them. 

The destructive scanner can be used to form the initial digital data set 

1 0 (IDDS) of the patient's teeth by milling and scanning a physical model, such as a plaster 
casting, of the teeth. To ensure a consistent orientation of the casting throughout the 
destructive scanning process, a scanning system operator pots the casting in potting 
material, such as Encase-It epoxy from CGI, and cures the material in a pressure vacuum 
(PV) chamber to form a mold. Placing the potting material into the PV chamber ensures 

1 5 that the material sets relatively quickly with virtually no trapped air bubbles. The color of 
the potting material is selected to contrast sharply with the color of the casting material to 
ensure the clarity of the scanned image. The operator then mounts the mold to a 
mounting plate and places the molding plate in the destructive scanning system. 

A slicing mechanism ("cutter") in the destructive scanning system mills a 

20 thin slice (typically between 0.001 " and 0.006" thick) from the mold, and a positioning 
arm places the milled surface near an optical scanner. The optical scanner, which may be 
an off-the-shelf device such as a flatbed scanner or a digital camera, scans the surface to 
create a 2D image data set representing the surface. The positioning arm then repositions 
the mold below the cutter, which again mills a thin slice from mold. The resulting output 

25 of the destructive scanning system is a 3D image data set, which later is converted into a 
digital model of surfaces, as described in detail below. A destructive scanning system 
and the corresponding destructive scanning and data processing are described in U.S. 
Patent 5,621,648. 

FIG. 12 shows a scanned image 1200 of an exposed surface of a plaster 
30 tooth casting 1202 embedded in an epoxy mold 1204. The black color of the epoxy mold 
1204 provides sharp contrast with the white color of the plaster casting 1202. An 
orientation guide 1206 appears in a corner of each image slice to ensure proper alignment 
of the image data after the destructive scan. The orientation guide 1206 can include a 
rigid structure, such as a piece of PVC tubing, embedded in the mold 1204. 




A patient's wax bite can be used to acquire the relative positions of the 
upper and lower teeth in centric occlusion. For a laser scan, this can be accomplished by 
first placing the lower cast in front of the scanner, with the teeth facing upwards, then 
placing the wax bite on top of the lower cast, and finally placing the upper cast on top of 
5 the lower cast, with the teeth facing downwards, resting on the wax bite. A cylindrical 
scan is then acquired for the lower and upper casts in their relative positions. The 
scanned data provides a digital model of medium resolution representing an object which 
is the combination of the patient's arches positioned in the same relative configuration as 
in the mouth. 

10 The digital model acts as a template guiding the placement of the two 

individual digital models (one per arch). More precisely, using software, for example the 
CyberWare alignment software, each digital arch is in turn aligned to the pair scan. The 
individual models are then positioned relative to each other corresponding to the arches in 
the patient's mouth. 

1 5 The waxbite can also be scanned separately to provide a second set of data 

about the teeth in the upper and lower arches. In particular, the plaster cast provides a 
"positive" image of the patient's teeth, from which one set of data is derived, and the 
waxbite provides a "negative" image of the teeth, from which a second, redundant set of 
data is derived. The two sets of data can then be matched to form a single data set 

20 describing the patient's teeth with increased accuracy and precision. The impression 
from which the plaster cast was made also can be used instead of or in addition to the 
waxbite. 

FIG. 13 is a flow chart of a process for deriving a single set of data from a 
positive data set and a negative data set. First, scan data representing positive and 

25 negative physical tooth models is obtained (steps 1300, 1302). If the scan data was 
acquired through a destructive scanning process, two digital 3D geometric models are 
constructed from the data, as described below (step 1304). Scan data acquired from an 
optical or range acquisition scanner system typically suffices as a geometric model. One 
of the geometric models is positioned to match roughly with the other model in the digital 

30 model space (step 1306), and an optimization process is performed to determine the best 
match between the models (step 1308). The optimization process attempts to match one 
point in each model to one point in the other model. Each pair of matched points then is 
combined into a single point to form a single set of data (step 1310). The combining of 
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matched points can be carried out in a variety of ways, for example, by averaging the 
coordinates of the points in each pair. 

While a laser scanning system typically must perform three scans to image 
a patient's full set of teeth adequately (one high resolution scan for each of the upper and 
5 lower casts and a medium resolution waxbite scan), the destructive scanning system 
described above can image a patient's full set of teeth adequately with only a single 
waxbite scan. Scanning both casts with the wax bite in place ensures that all important 
surfaces of the upper and lower teeth are captured during a destructive scan. Scanning 
both casts in this manner also provides a high resolution image data set that preserves the 

1 0 relation between the patient's upper and lower teeth. Like the potting material described 
above, the wax bite should have a color that contrasts sharply with the color of the casting 
material to ensure the clarity of the scanned image. The wax bite may be the same color 
as the potting material if no contrast is desired between the waxbite and the potting 
material. Alternatively, the color of the wax bite may contrast sharply with the tooth 

1 5 casts and the potting material if an image of the wax bite is needed. 

In addition to the 3D image data gathered by laser scanning or destructive 
scanning the exposed surfaces of the teeth, a user may wish to gather data about hidden 
features, such as the roots of the patient's teeth and the patient's jaw bones. This 
information is used to build a more complete model of the patient's dentition and to show 

20 with more accuracy and precision how the teeth will respond to treatment. For example, 
information about the roots allows modeling of all tooth surfaces, instead of just the 
crowns, which in turn allows simulation of the relationships between the crowns and the 
roots as they move during treatment. Information about the patient's jaws and gums also 
enables a more accurate model of tooth movement during treatment. For example, an x- 

25 ray of the patient's jaw bones can assist in identifying ankylose teeth, and an MRI can 

provide information about the density of the patient's gum tissue. Moreover, information 
about the relationship between the patient's teeth and other cranial features allows 
accurate alignment of the teeth with respect to the rest of the head at each of the treatment 
steps. Data about these hidden features may be gathered from many sources, including 

30 2D and 3D x-ray systems, CT scanners, and magnetic resonance imaging (MRI) systems. 
Using this data to introduce visually hidden features to the tooth model is described in 
more detail below. 

Developing an orthodontic treatment plan for a patient involves 
manipulating the IDDS at a computer or workstation having a suitable graphical user 
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interface (GUI) and software appropriate for viewing and modifying the images. Specific 
aspects of the software will be described in detail hereinafter. However, dental 
appliances having incrementally differing geometries can be produced by non-computer- 
aided techniques. For example, plaster casts obtained as described above may be cut 
5 using knives, saws, or other cutting tools in order to permit repositioning of individual 
teeth within the casting. The disconnected teeth may then be held in place by soft wax or 
other malleable material, and a plurality of intermediate tooth arrangements can then be 
prepared using such a modified plaster casting of the patient's teeth. The different 
arrangements can be used to prepare sets of multiple appliances, generally as described 
1 0 below, using pressure and vacuum molding techniques. 

Creating a 3D surface model of the teeth 

Many types of scan data, such as that acquired by an optical scanning 
system, provide a 3D geometric model (e.g., a triangular surface mesh) of the teeth when 
acquired. Other scanning techniques, such as the destructive scanning technique 

1 5 described above, provide data in the form of volume elements ("voxels") that can be 
converted into a digital geometric model of the tooth surfaces. 

FIG. 14 is a flowchart of a process for forming a surface mesh from voxel 
image data. This approach involves receiving the image data from the destructive scanner 
(step 1400), processing the data to isolate the object to be modeled (step 1401), and 

20 applying a conventional "marching cubes" technique to create a surface mesh of the 
object (step 1402). 

Each set of image data can include images of multiple tooth casts or of a 
tooth cast and extraneous, "noisy" objects, such as air bubbles in the potting material. 
The system identifies each object in the image by assigning each voxel a single-bit binary 

25 value (e.g., "0" for black and "1" for white) based on the voxel's 8-bit gray scale image 
value, and then connecting adjacent voxels that have been assigned the same single-bit 
value. Each group of connected voxels represents one of the objects in the image. The 
system then isolates the tooth casting to be modeled by masking from the image all 
objects except the tooth casting of interest. The system removes noise from the masked 

30 image data by passing the data through a low-pass filter. 

Once the tooth casting is isolated from the image data, the system 
performs a conventional marching cubes technique to locate tooth and tissue surfaces in 
the image data. This technique involves the identification of pairs of adjacent voxels 
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having 8-bit image values that fall on opposite sides of a selected threshold value. 
Specifically, each voxel has an associated image value, typically a number between 0 and 
255, that represents the image color or gray scale value at that voxel. Because the tooth 
casting and the surrounding potting material have sharply contrasting colors (see FIG. 
5 12), the image values at voxels forming the image of the tooth casting differ greatly from 
the values at voxels forming the image of the surrounding material. Therefore, the 
marching cube algorithm can locate the tooth surfaces by identifying the voxels at which 
a sharp transition in image value occurs. The algorithm can position the surface precisely 
between two voxels by determining the difference between the threshold value and the 

1 0 image value at each voxel and then placing the surface a corresponding distance from the 
centerpoint of each voxel. 

In some implementations, after the marching cubes algorithm is applied, 
the resulting mesh undergoes a smoothing operation to reduce the jaggedness on the 
surfaces of the tooth model caused by the marching cubes conversion (step 1404). A 

15 conventional smoothing operation can be used, such as one that moves individual triangle 
vertices to positions representing the averages of connected neighborhood vertices to 
reduce the angles between triangles in the mesh. 

Another optional step is the application of a decimation operation to the 
smoothed mesh to eliminate data points, which improves processing speed (step 1406). 

20 Conventional decimation operations identify pairs of triangles that lie almost in the same 
plane and combine each identified pair into a single triangle by eliminating the common 
vertex. The decimation operation used here also incorporates orthodontic-specific 
decimation rules, which rely on an understanding of the general characteristics of the 
teeth and gums and of the orthodontic appliances that will be used to carry out a treatment 

25 plan. For example, aligners typically do not contact the portions of the tooth surfaces 

adjacent to the gums, so these tooth surfaces can be modeled less accurately than the rest 
of the tooth. The decimation operation incorporates this knowledge by decimating more 
heavily along the gum line. When an appliance such as a polymeric shell aligner will be 
used to treat the patient's teeth, the algorithm also decimates more heavily on the sides of 

30 the teeth, where the aligner usually is required only to push orthogonally to the surface, 
than it decimates on the tops of the teeth, where the aligner usually must form a solid 

grip- 
After the smoothing and decimation operation are performed, an error 
value is calculated based on the differences between the resulting mesh and the original 
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mesh or the original data (step 1408), and the error is compared to an acceptable threshold 
value (step 1410). The smoothing and decimation operations are applied to the mesh 
once again if the error does not exceed the acceptable value. The last set of mesh data 
that satisfies the threshold is stored as the tooth model (step 1412). 

5 Creating 3D models of the individual teeth 

Once a 3D model of the tooth surfaces has been constructed, models of the 
patient's individual teeth can be derived. In one approach, individual teeth and other 
components are "cut" to permit individual repositioning or removal of teeth in or from the 
digital data. After the components are "freed," a prescription or other written 
10 specification provided by the treating professional is followed to reposition the teeth. 
Alternatively, the teeth may be repositioned based on the visual appearance or based on 
rules and algorithms programmed into the computer. Once an acceptable final 
arrangement has been created, the final tooth arrangement is incorporated into a final 
digital data set (FDDS). 

15 Based on both the IDDS and the FDDS, a plurality of intermediate digital 

data sets (INTDDS's) are generated to correspond to successive intermediate tooth 
arrangements. The system of incremental position adjustment appliances can then be 
fabricated based on the INTDDS's, as described in more detail below. Segmentation of 
individual teeth from the tooth model and determining the intermediate and final positions 

20 of teeth also are described in more detail below. 

Simplifying the 3D model 

FIG. 3 illustrates a representative technique for user-assisted manipulation 
of the IDDS to produce the FDDS on the computer. Usually, the data from the digital 
scanner is acquired in high resolution form. In order to reduce the computer time 

25 necessary to generate images, a parallel set of digital data representing the IDDS at a 
lower resolution can be created. The user can manipulate the lower resolution images 
while the computer updates the high resolution data set as necessary. The user can also 
view and manipulate the high resolution model if the extra detail provided in that model 
is useful. The IDDS also can be converted into a quad edge data structure if not already 

30 present in that form. A quad edge data structure is a standard topological data structure 
defined in Primitives for the Manipulation of General Subdivisions and the Computation 
ofVoronoi Diagrams, ACM Transactions of Graphics, Vol. 4, No. 2, April 1985, pp. 74- 
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123. Other topological data structures, such as the winged-edge data structure, could also 
be used. 

As an initial step, while viewing the three-dimensional image of the 
patient's jaw. including the teeth, gingivae, and other oral tissue, the user usually deletes 
5 structure which is unnecessary for image manipulation and final production of an 

appliance. These unwanted sections of the model may be removed using an eraser tool to 
perform a solid modeling subtraction. The tool is represented by a graphic box. The 
volume to be erased (the dimensions, position, and orientation of the box) are set by the 
user employing the GUI. Typically, unwanted sections would include extraneous gum 

1 0 area and the base of the originally scanned cast. Another application for this tool is to 
stimulate the extraction of teeth and the "shaving down" of tooth surfaces. This is 
necessary when additional space is needed in the jaw for the final positioning of a tooth to 
be moved. The treating professional may choose to determine which teeth will be shaved 
and which teeth will be extracted. Shaving allows the patient to maintain teeth when only 

1 5 a small amount of space is needed. Typically, extraction and shaving are used in the 

treatment planning only when the actual patient teeth are to be extracted or shaved prior 
to initiating repositioning. 

Removing unwanted or unnecessary sections of the model increases data 
processing speed and enhances the visual display. Unnecessary sections include those not 

20 needed for creation of the tooth repositioning appliance. The removal of these unwanted 
sections reduces the complexity and size of the digital data set, thus accelerating 
manipulations of the data set and other operations. 

After the user positions and sizes the eraser tool and instructs the software 
to erase the unwanted section, all triangles within the box set by the user are removed and 

25 the border triangles are modified to leave a smooth, linear border. The software deletes 
all of the triangles within the box and clips all triangles which cross the border of the box. 
This requires generating new vertices on the border of the box. The holes created in the 
model at the faces of the box are retriangulated and closed using the newly created 
vertices. 

30 In alternative embodiments, the computer automatically simplifies the 

digital model by performing the user-oriented functions described above. The computer 
applies a knowledge of orthodontic relevance to determine which portions of the digital 
model are unnecessary for image manipulation. 
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SEGMENTING THE TEETH IN THE 3D MODEL 
Human-assisted segmentation 

The saw tool is used to define the individual teeth (or possibly groups of 
teeth) to be moved. The tool separates the scanned image into individual graphic 
5 components enabling the software to move the tooth or other component images 

independent of remaining portions of the model. In one embodiment, the saw tool defines 
a path for cutting the graphic image by using two cubic B-spline curves lying in space, 
possibly constrained to parallel planes, either open or closed. A set of lines connects the 
two curves and shows the user the general cutting path. The user may edit the control 

10 points on the cubic B-splines, the thickness of the saw cut, and the number of erasers 
used, as described below. 

In an alternative embodiment, the teeth are separated by using the saw as a 
"coring" device, cutting the tooth from above with vertical saw cuts. The crown of the 
tooth, as well as the gingivae tissue immediately below the crown are separated from the 

1 5 rest of the geometry, and treated as an individual unit, referred to as a tooth. When this 
model is moved, the gingivae tissue moves relative to the crown, creating a first order 
approximation of the way that the gingivae will reform within a patient's mouth. 

Each tooth may also be separated from the original trimmed model. 
Additionally, a base may be created from the original trimmed model by cutting off the 

20 crowns of the teeth. The resulting model is used as a base for moving the teeth. This 
facilitates the eventual manufacture of a physical mold from the geometric model, as 
described below. 

Thickness: When a cut is used to separate a tooth, the user will usually 
want the cut to be as thin as possible. However, the user may want to make a thicker cut, 

25 for example, when shaving down surrounding teeth, as described above. Graphically, the 
cut appears as a curve bounded by the thickness of the cut on one side of the curve. 

Number of Erasers: A cut is comprised of multiple eraser boxes arranged 
next to each other as a piecewise linear approximation of the Saw Tool's curve path. The 
user chooses the number of erasers, which determines the sophistication of the curve 

30 created: the greater the number of segments, the more accurately the cutting will follow 
the curve. The number of erasers is shown graphically by the number of parallel lines 
connecting the two cubic B-spline curves. Once a saw cut has been completely specified 
the user applies the cut to the model. The cut is performed as a sequence of erasings, as 
shown in FIG. 4A. FIG. 4B shows a single erasing iteration of the cut as described in the 
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algorithm for a open ended B-spline curve. For a vertical cut. the curves are closed, with 
P A [O] and P A [S] being the same point and P B [O] and P B [S] being the same point. 

In one embodiment the software may automatically partition the saw tool 
into a set of erasers based upon a smoothness measure input by the user. The saw is 
5 adaptively subdivided until an error metric measures the deviation from the ideal 

representation to the approximate representation to be less than a threshold specified by 
the smoothness setting. One error metric compares the linear length of the subdivided 
curve to the arclength of the ideal spline curve. When the difference is greater than a 
threshold computed from the smoothness setting, a subdivision point is added along the 
1 0 spline curve. 

A preview feature may also be provided in the software. The preview 
feature visually displays a saw cut as the two surfaces that represent opposed sides of the 
cut. This allows the user to consider the final cut before applying it to the model data set. 

After the user has completed all desired cutting operations with the saw 
15 tool, multiple graphic solids exist. However, at this point, the software has not 

determined which triangles of the quad edge data structure belong to which components. 
The software chooses a random starting point in the data structure and traverses the data 
structure using adjacency information to find all of the triangles that are attached to each 
other, identifying an individual component. This process is repeated starting with the 
20 triangle whose component is not yet determined. Once the entire data structure is 
traversed, all components have been identified. 

To the user, all changes made to the high resolution model appear to occur 
simultaneously in the low resolution model, and vice versa. However, there is not a one- 
to-one correlation between the different resolution models. Therefore, the computer 
25 "matches" the high resolution and low resolution components as best as it can subject to 
defined limits. One process for doing so is described in FIG. 5. 

Automated segmentation 

The system can optionally include a segmentation subsystem that performs 
automatic or semi-automatic segmentation of the 3D dentition model into models of 
30 individual teeth. The segmentation subsystem is advantageously implemented as one or 
more computer program processes implementing a segmentation process. In alternative 
implementations, the segmentation process can act on the 3D volume data or on the 3D 
surface mesh. The segmentation process applies conventional feature detection 
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techniques tailored to exploit the characteristics and known features of teeth. For 
example, feature detection algorithms generally act on images in which the features to be 
distinguished from each other have different colors or shades of gray. Features to be 
detected also usually are separated spatially from each other. However, features to be 
5 detected in a 2D or 3D image of a plaster tooth casting (e.g. ? the individual teeth and the 
gum tissue) all have the same color (white), and some features, such as an individual 
tooth and the surrounding gum tissue, have no spacial separation. 

The segmentation process can be implemented to employ any of several 
feature detection techniques and advantageously uses a combination of techniques to 

10 increase the accuracy of feature identification. One feature detection technique uses color 
analysis to distinguish objects based on variations in color. Color analysis can be used in 
situations where individual teeth are separated by gaps large enough for the potting 
material to fill. Because the tooth casting and the potting material have contrasting 
colors, these teeth appear in the model as white areas separated by thin strips of black. 

15 Another feature detection technique uses shape analysis to distinguish 

certain features, such as tooth from gum. In general, tooth surfaces are smooth while gum 
surfaces have texture, and the teeth and gums typically form a U-shaped ridge where they 
meet. Detecting these features through shape analysis assists in distinguishing tooth from 
gum. Shape analysis can also detect individual teeth, for example by searching for the 

20 largest objects in the 3D image or by recognizing the cusps of a molar as four isolated 
patches of one color arranged in a certain pattern. One cusp-detection algorithm is 
described below. 

Other feature detection techniques use databases of known cases or 
statistical information against which a particular 3D image is matched using conventional 

25 image pattern matching and data fitting techniques. One such technique, known as 

"Maximum a posteriori" (MAP), uses prior images to model pixel values corresponding 
to distinct object types (classes) as independent random variables with normal (Gaussian) 
distributions whose parameters (mean and variance) are selected empirically. For each 
class, a histogram profile is created based on a Gaussian distribution with the specified 

30 mean and variance. The prior images supply for each pixel and each class the probability 
that the pixel belongs to the class, a measure which reflects the relative frequency of each 
class. Applying Bayes' Rule to each class, the pixel values in the input image are scaled 
according to the prior probabilities, then by the distribution function. The result is a 
posterior probability that each pixel belongs to each class. The Maximum a posteriori 
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(MAP) approach then selects for each pixel the class with the highest posterior 
probability as the output of the segmentation. 

Another feature detection technique uses automatic detection of tooth 
cusps. Cusps are pointed projections on the chewing surface of a tooth. In one 
5 implementation cusp detection is performed in two stages: (1) a "detection" stage, during 
which a set of points on the tooth are determined as candidates for cusp locations; and (2) 
a "rejection" stage, during which candidates from the set of points are rejected if they do 
not satisfy a set of criteria associated with cusps. 

One process for the "detection" stage is set forth in FIG. 6A. In the 
10 detection stage, a possible cusp is viewed as an "island" on the surface of the tooth, with 
the candidate cusp at the highest point on the island. "Highest" is measured with respect 
to the coordinate system of the model, but could just as easily be measured with respect to 
the local coordinate system of each tooth if detection is performed after the cutting phase 
of treatment. 

1 5 The set of all possible cusps is determined by looking for all local maxima 

on the tooth model that are within a specified distance of the top of the bounding box of 
the model. First, the highest point on the model is designated as the first candidate cusp. 
A plane is passed through this point, perpendicular to the direction along which the height 
of a point is measured. The plane is then lowered by a small predetermined distance 

20 along the Z axis. Next, all vertices connected to the tooth and which are above the plane 
and on some connected component are associated with the candidate cusp as cusps. This 
step is also referred to as the "flood fill" step. From each candidate cusp point, outward 
"flooding" is performed, marking each vertex on the model visited in this matter as "part 
of the corresponding candidate cusp. After the flood fill step is complete, every vertex 

25 on the model is examined. Any vertex that is above the plane and has not been visited by 
one of the flood fills is added to the list of candidate cusps. These steps are repeated until 
the plane is traveled a specified distance. 

While this iterative approach can be more time consuming than a local 
maximum search, the approach described above leads to a shorter list of candidate cusps. 

30 Since the plane is lowered a finite distance at each step, very small local maxima that can 
occur due to noisy data are skipped over. 

After the "detection" stage, the cusp detection algorithm proceeds with the 
"rejection" stage. One process for the "rejection" stage is set forth in FIG. 6B. In this 
stage, the local geometries around each of cusp candidates are analyzed to determine if 
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they possess M non-cusp-like features." Cusp candidates that exhibit "non-cusp-like 
features" are removed from the list of cusp candidates. 

Various criteria may be used to identify "non-cusp-like features." 
According to one test, the local curvature of the surface around the cusp candidate is used 
5 to determine whether the candidate possesses non-cusp-like features. As depicted in FIG. 
6B, the local curvature of the surface around the cusp candidate is approximated and then 
analyzed to determine if it is too large (very pointy surface) or too small (very flat 
surface), in which case the candidate is removed from the list of cusp candidates. 
Conservative values are used for the minimum and maximum curvature values to ensure 

1 0 that genuine cusps are not rejected by mistake. 

Under an alternative test, a measure of smoothness is computed based on 
the average normal in an area around the candidate cusp. If the average normal deviates 
from the normal at the cusp by more than a specified amount, the candidate cusp is 
rejected. In one embodiment, the deviation of a normal vector N from the cusp normal 

15 CN is approximated by the formula: 

deviation = 1 - Abs(N-CN), 
which is zero at no deviation, and 1 when N and CN are perpendicular. 

For both the human-assisted and automated segmentation techniques, the 
clinician can simplify the tooth identification process by marking the physical tooth 

20 model before the model is scanned. Upon scanning, these marks become part of the 
digital tooth model. The types of marks that the clinician might use include marks 
identifying the rotational axis of a tooth, marks identifying the principal axis of a tooth 
(e.g., a straight line marked on the tooth's occlusal edge), and marks identifying the 
boundaries between teeth. A mark identifying the rotational axis of a tooth often is used 

25 to restrict how the tooth can rotate during the course of treatment. The clinician also may 
wish to paint the teeth in the physical model with various colors to assist in segmenting 
the individual teeth from the digital tooth model. 

Adding roots and hidden tooth surfaces to the individual tooth models 

The system can optionally be configured to add roots and hidden surfaces 
30 to the tooth models to allow more thorough and accurate simulation of tooth movement 
during treatment. In alternative implementations, this information is added automatically 
without human assistance, semi-automatically with human assistance, or manually by 
human operator, using a variety of data sources. 
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In some implementations, 2D and 3D imaging systems, such as x-ray 
systems, computed tomography (CT) scanners, and MRI systems, are used to gather 
information about the roots of the patient's teeth. For example, several 2D x-ray images 
of a tooth taken in different planes allow the construction of a 3D model of the tooth's 
5 roots. Information about the roots is available by visual inspection of the x-ray image and 
by application of a computer-implemented feature identification algorithm to the x-ray 
data. The system adds the roots to the tooth model by creating a surface mesh 
representing the roots. Physical landmarks on the patient's teeth, e.g., cavities or cusps, 
are extracted from the 2D and 3D data and are used to register the roots to the tooth 

10 model. Like the roots, these landmarks can be extracted manually or by use of a feature 
detection algorithm. 

Another alternative for the addition of roots and hidden surfaces is to 
model typical root and crown shapes and to modify the digital model of each tooth to 
include a root or a hidden surface corresponding to a typical shape. This approach 

1 5 assumes that the roots and hidden surfaces of each patient's teeth have typical shapes. A 
geometric model of each typical shape is acquired, e.g., by accessing an electronic 
database of typical root and crown models created before the analysis of a particular 
patient's teeth begins. Portions of the typical root and crown models are added to the 
individual tooth models as needed to complete the individual tooth models. 

20 Yet another alternative for the addition of roots and hidden surfaces is the 

extrapolation of the 3D tooth model to include these features based on observed 
characteristics of the tooth surfaces. For example, the system can use the curvature of a 
particular molar between the tips of the cusps and the gumline to predict the shape of the 
roots for that molar. In other implementations, x-ray and CT scan data of a patient's teeth 

25 are used to provide comparison points for extrapolation of the patients roots and hidden 
surfaces. Models of typical root and crown shapes also can be used to provide 
comparison points for root and hidden surface extrapolation. 

DETERMINING THE FINAL TOOTH POSITIONS 

Once the teeth have been separated, the FDDS can be created from the 
30 IDDS. The FDDS is created by following the orthodontists' prescription to move the 

teeth in the model to their final positions. In one embodiment, the prescription is entered 
into a computer, which automatically computes the final positions of the teeth. In 
alternative other embodiments, a user moves the teeth into their final positions by 
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independently manipulating one or more teeth while satisfying the constraints of the 
prescription. Various combinations of the above described techniques may also be used 
to arrive at the final tooth positions. 

One method for creating the FDDS involves moving the teeth in a 
5 specified sequence. First, the centers of each of the teeth are aligned to a standard arch. 
Then, the teeth are rotated until their roots are in the proper vertical position. Next, the 
teeth are rotated around their vertical axis into the proper orientation. The teeth are then 
observed from the side, and translated vertically into their proper vertical position. The 
inclusion of roots in the tooth models, described more fully above, ensures' the proper 

10 vertical orientation of the entire tooth, not just the crown. Finally, the two arches are 
placed together, and the teeth moved slightly to ensure that the upper and lower arches 
properly mesh together. The meshing of the upper and lower arches together is visualized 
using a collision detection process to highlight the contacting points of the teeth in red. 

Apart from its role in identifying individual teeth, cusp detection also is 

15 useful in determining final tooth orientation. For example, the cusps on a typical molar 
are relatively level when the tooth is oriented vertically, so the relative positions of the 
cusp tips indicate the tooth's position. Cusp detection therefore is incorporated into the 
final position determination. 

One tool for use in visualizing the interaction of a patient's upper and 

20 lower teeth at the final positions is a computer-implemented "virtual" articulator. The 

virtual articulator provides a graphical display that simulates the operation of the patient's 
jaw or the operation of a conventional mechanical articulator attached to a physical model 
of the patient's teeth. In particular, the virtual articulator orients the digital models of the 
patient's upper and lower arches in the same manner that the patient's physical arches 

25 will be oriented in the patient's mouth at the end of treatment. The articular then moves 
the arch models through a range of motions that simulate common motions of the human 
jaw. 

The quality of the virtual articulator's simulation depends on the types of 
information used to create the articulator and the tooth models. In some implementations, 
30 the virtual articulator includes a digital model of a mechanical articular created, for 
example, from a computer-aided drafting (CAD) file or image data gathered during a 
laser scan of the articulator. Other implementations include a digital model of human 
jaws created, for example, from 2D or 3D x-ray data, CT scan data, or mechanical 
measurements of the jaws, or from a combination of these types of data. In many 
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respects, the most useful virtual articulator is the one that simulates the jaws of the patient 
whose teeth are being treated, which is created from image data or mechanical 
measurements of the patient's head. 

Animation instructions define the movements that the virtual articulator 
5 simulates. Like the articulator itself, the animation instructions are derived from a variety 
of sources. The animation instructions associated with the simulation of a mechanical 
articulator require little more than a mathematical description of the motion of a 
mechanical hinge. A virtual articulator simulating the human jaw, on the other hand, 
requires a more complex set of instructions, based on human anatomical data. One 

10 method of building this set of instructions is the derivation of mathematical equations 
describing the common motions of an ideal human jaw. Another method is through the 
use of a commercially available jaw-tracking system, which contacts a person's face and 
provides digital information describing the motion of the lower jaw. X-ray and CT scan 
data also provide information indicating how the teeth and jaws relate to each other and to 

15 the rest of the person's head. Jaw-tracking systems and x-ray and CT scan data are 
particularly useful in developing an articulator that simulates a particular patient's 
anatomy. 

As the virtual articulator simulates the motion of a patient's teeth, a 
collision detection process, such as one implementing the oriented bounding box (OBB) 
20 algorithm described below, determines whether and how the patient's teeth will collide 

during the normal course of oral motion. Visual indicators, such as red highlights, appear 
on a displayed image of the teeth to indicate collision points. The final tooth positions are 
adjusted, automatically or manually, to avoid collisions detected by the collision detection 
algorithm. 

25 An automated system for determining final tooth positions and creating the 

FDDS is described in the above-mentioned U.S. patent application 

(attorney docket number 09943/003001). That application describes a computer- 
implemented process for generating a set of final positions for a patient's teeth. The 
process involves creating an ideal model of final tooth positions based on "ideal" tooth 

30 arrangements, repositioning the individual teeth in a digital model of the patient's teeth to 
mimic the ideal model, and modeling the motion of the patient's jaw to perfect the final 
tooth arrangement. 

The display and use of orthodontic-specific information also assists in the 
determination of final tooth positions. In some implementations, a user can elect to have 
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malocclusion indices, such as Peer Assessment Review (PAR) metrics, shape-based 
metrics, or distance-based metrics, calculated and displayed with the final tooth positions. 
If the user is not satisfied with values of the displayed index or indices, the user can adjust 
the final tooth positions manually until the parameters fall within acceptable ranges. If 
5 the tooth positioning system is fully automated, the orthodontic-specific parameters are 
provided as feedback and used to adjust the final tooth positions until the parameters fall 
within acceptable ranges. 

For human-assisted tooth positioning, the human user also can elect to 
have positioning tips displayed. Tips are available, for example, to suggest a treatment 
10 path for an individual tooth and to warn of excessive forces that might cause patient 
discomfort or compromise the mechanical integrity of the orthodontic appliance. Tips 
also are available to suggest target positions that best suit the patient's jaw structure and 
that ensure proper inner digitation and occlusion parameters. 

DETERMINING THE STEPS FROM THE INITIAL TO THE FINAL POSITION 
1 5 After the teeth and other components have been placed or removed to 

produce a model of the final tooth arrangement, it is necessary to generate a treatment 
plan, as illustrated in FIG. 7. The treatment plan will ultimately produce the series of 
INTDDS's and FDDS as described previously. To produce these data sets, it is necessary 
to define or map the movement of selected individual teeth from the initial position to the 
20 final position over a series of successive steps. In addition, it may be necessary to add 
other features to the data sets in order to produce desired features in the treatment 
appliances. For example, it may be desirable to add wax patches to the image in order to 
define cavities or recesses for particular purposes, such as to maintain a space between 
the appliance and particular regions of the teeth or jaw in order to reduce soreness of the 
25 gums, avoid periodontal problems, allow for a cap, and the like. Additionally, it will 

often be necessary to provide a receptacle or aperture intended to accommodate an anchor 
which is to be placed on a tooth in order to permit the tooth to be manipulated in a 
manner that requires the anchor, e.g., to be lifted relative to the jaw. 

Accounting for physical constraints and additions to the model 
30 Some methods for manufacturing the tooth repositioning appliances 

require that the separate, repositioned teeth and other components be unified into a single 
continuous structure in order to permit manufacturing. In these instances, "wax patches 
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are used to attach otherwise disconnected components of the INTDDS's. These patches 
are added to the data set underneath the teeth and above the gum so that they do not effect 
the geometry of the tooth repositioning appliances. The application software provides for 
a variety of wax patches to be added to the model, including boxes and spheres with 
5 adjustable dimensions. The wax patches that are added are treated by the software as 
additional pieces of geometry, identical to all other geometries. Thus, the wax patches 
can be repositioned during the treatment path, as can the teeth and other components. 
One method of separating the teeth using vertical coring, as described above, removes the 
need for most of these "wax patches". 

10 In the manufacturing process, adding a wax patch to the graphic model 

will generate a positive mold that has the same added wax patch geometry. Because the 
mold is a positive of the teeth and a polymeric appliance is a negative of the teeth, when 
the appliance is formed over the mold, the appliance will also form around the wax patch 
that has been added to the mold. When placed in the patient's mouth, the appliance will 

1 5 thus allow for a space between the inner cavity surface of the appliance and the patient's 
teeth or gums. Additionally, the wax patch may be used to form a recess or aperture 
within the appliance which engages an anchor placed on the teeth in order to move the 
tooth in directions which could not otherwise be accomplished. 

For some patients an optimal treatment plan requires the interaction of 

20 aligners with tooth attachments, such as brackets and anchors, to ensure proper 

orthodontic correction in a reasonable amount of time. In these situations, the aligners 
must grip the attachments to ensure that the proper force is exerted on the patient's teeth. 
For example, an aligner may be designed to grip an anchor planted in the patient's jaw to 
move the patient's teeth back in the jaw. Likewise, an aligner may grip a bracket attached 

25 to a patient's tooth to increase the aligner's leverage or grip on the tooth. 

The creation of digital attachment models allows the system to model the 
effects of attachments in analyzing the digital model of a patient's teeth. Each attachment 
model represents a physical attachment that may be placed in a patient's mouth, generally 
on a tooth, during the course of treatment. Many attachments, such as conventional 

30 brackets, are available in standard shapes and sizes, the models of which can be selected 
from a library of virtual attachments and added to a patient's tooth model. Other 
attachments are patient-specific and must be modeled by the user for inclusion in the 
digital tooth model. The presence of virtual attachments in a patient's tooth model 
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ensures that the aligners fabricated for the patient's treatment plan will accommodate the 
corresponding physical attachments placed in the patient's mouth during treatment. 

The wax patches and virtual attachments described above, and individual 
components of the tooth model, can be reduced or enlarged in size, which will result in a 
5 manufactured appliance having a tighter or looser fit. 

Selecting the intermediate treatment stages . 

Number of Treatment Stages: The user can change the number of desired 
treatment stages from the initial to the target states of the teeth. Any component that is 
not moved is assumed to remain stationary, and thus its final position is assumed to be the 

10 same as the initial position (likewise for all intermediate positions, unless one or more 
key frames are defined for that component). 

Key frames: The user may also specify "key frames" by selecting an 
intermediate state and making changes to component position(s). In some embodiments, 
unless instructed otherwise, the software automatically linearly interpolates between all 

1 5 user-specified positions (including the initial position, all key frame positions, and the ' 
target position). For example, if only a final position is defined for a particular 
component, each subsequent stage after the initial stage will simply show the component 
an equal linear distance and rotation (specified by a quaternion) closer to the final 
position. If the user specifies two key frames for that component, the component will 

20 "move" linearly from the initial position through different stages to the position defined 
by the first key frame. It will then move, possibly in a different direction, linearly to the 
position defined by the second key frame. Finally, it will move, possibly in yet a 
different direction, linearly to the target position. 

These operations may be done independently to each component, so that a 

25 key frame for one component will not affect another component, unless the other 

component is also moved by the user in that key frame. One component may accelerate 
along a curve between one pair of stages (e.g., stages 3 and 8 in a treatment plan having 
that many stages), while another moves linearly between another pair of stages (e.g., 
stages 1 to 5), and then changes direction suddenly and slows down along a linear path to 

30 a later stage (e.g., stage 10). This flexibility allows a great deal of freedom in planning a 
patient's treatment. 

In some implementations, non-linear interpolation is used instead of or in 
addition to linear interpolation to construct a treatment path among key frames. In 
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general, a non-linear path, such as a spline curve, created to fit among selected points is 
shorter than a path formed from straight line segments connecting the points. A 
"treatment path" describes the transformation curve applied to a particular tooth to move 
the tooth from its initial position to its final position. A typical treatment path includes 
5 some combination of rotational and translational movement of the corresponding tooth, as 
described above. 

FIGS. 15A and 15B show a linearly interpolated treatment path and a non- 
linearly interpolated path, respectively, connecting an initial tooth position I to a final 
tooth position F. The linearly interpolated path consists of straight line segments 

1 0 connecting the initial and final tooth positions as well as the four intermediate tooth 
positions Ii, h, The non-linear interpolated path consists of a curved line fitted 
among the intermediate tooth positions. The curved path can be formed using a 
conventional spline-curve fitting algorithm. 

FIG. 16 is a flow chart of a computer-implemented process for generating 

1 5 non-linear treatment paths along which a patient's teeth will travel during treatment. The 
non-linear paths usually are generated automatically by computer program, in some cases 
with human assistance. The program receives as input the initial and final positions of the 
patient's teeth and uses this information to select intermediate positions for each tooth to 
be moved (step 1600). The program then applies a conventional spline curve calculation 

20 algorithm to create a spline curve connecting each tooth's initial position to the tooth's 
final position (step 1602). In many situations, the curve is constrained to follow the 
shortest path between the intermediate positions. The program then samples each spline 
curve between the intermediate positions (step 1604) and applies the collision detection 
algorithm to the samples (step 1606). If any collisions are detected, the program alters 

25 the path of at least one tooth in each colliding pair by selecting a new position for one of 
. the intermediate steps (step 1608) and creating a new spline curve (1602). The program 
then samples the new path (1604) and again applies the collision detection algorithm 
(1606). The program continues in this manner until no collisions are detected. The 
routine then stores the paths, e.g., by saving the coordinates of each point in the tooth at 

30 each position on the path in an electronic storage device, such as a hard disk (step 1610). 

The path-generating program, whether using linear or non-linear 
interpolation, selects the treatment positions so that the tooth's treatment path has 
approximately equal lengths between each adjacent pair of treatment steps. The program 
also avoids treatment positions that force portions of a tooth to move with more than a 
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given maximum velocity. FIG. 15C shows a tooth that is scheduled to move along a first 
path Tl from an initial position Tl i to a final position TI3 through an intermediate 
position Tl 2 , which lies closer to the final position Tl 3 . Another tooth is scheduled to 
move along a shorter path T2 from an initial position T2\ to a final position T2 3 through 
5 an intermediate position T2 2? which is equidistant from the initial and final positions T2i, 
T23. In this situation, the program may choose to insert a second intermediate position 
Tl 4 along the first path Tl that is approximately equidistant from the initial position Tl 1 
and the intermediate position Tl 2 and that is separated from these two positions by 
approximately the same distance that separates the intermediate position Tl 2 from the 

10 final position Tl 3 . 

Altering the first path Tl in this manner ensures that the first tooth will 
move in steps of equal size. However, altering the first path Tl also introduces an 
additional treatment step having no counterpart in the second path T2. The program can 
respond to this situation in a variety of ways, such as by allowing the second tooth to 

1 5 remain stationary during the second treatment step (i.e., as the first tooth moves from one 
intermediate position Tl 4 to the other intermediate position TI3) or by altering the second 
path T2 to include four equidistant treatment positions. The program determines how to 
respond by applying a set of orthodontic constraints , that restrict the movement of the 
teeth. 

20 Orthodontic constraints that may be applied by the path-generating 

program include the minimum and maximum distances allowed between adjacent teeth at 
any given time, the maximum linear or rotational velocity at which a tooth should move, 
the maximum distance over which a tooth should move between treatment steps, the 
shapes of the teeth, the characteristics of the tissue and bone surrounding the teeth (e.g., 

25 ankylose teeth cannot move at all), and the characteristics of the aligner material (e.g., the 
maximum distance that the aligner can move a given tooth over a given period of time). 
For example, the patient's age and jaw bone density may dictate certain "safe limits" 
beyond which the patient's teeth should not forced to move. In general, a gap between 
two adjacent, relatively vertical and non-tipped central and lateral teeth should not close 

30 by more than about 1 mm every seven weeks. The material properties of the orthodontic 
appliance also limit the amount by which the appliance can move a tooth. For example, 
conventional retainer materials usually limit individual tooth movement to approximately 
0.5 mm between treatment steps. The constraints have default values that apply unless 
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patient-specific values are calculated or provided by a user. Constraint information is 
available from a variety of sources, including text books and treating clinicians. 

In selecting the intermediate positions for each tooth, the path-generating 
program invokes the collision detection program to determine whether collisions will 
5 occur along the chosen paths. The program also inspects the patient's occlusion at each 
treatment step along the path to ensure that the teeth align to form an acceptable bite 
throughout the course of treatment. If collisions or an unacceptable bite will occur, or if a 
required constraint cannot be satisfied, the program iteratively alters the offending tooth 
path until all conditions are met. The virtual articulator described above is one tool for 

1 0 testing bite occlusion of the intermediate treatment positions. 

As shown in FIG. 17, once the path-generating program has established 
collision-free paths for each tooth to be moved, the program calls an optimization routine 
that attempts to make the transformation curve for each tooth between the initial and final 
positions more linear. The routine begins by sampling each treatment path at points 

15 between treatment steps (step 1702), e.g., by placing two sample points between each 

treatment step, and calculating for each tooth a more linear treatment path that fits among 
the sample points (step 1704). The routine then applies the collision detection algorithm 
to determine whether collisions result from the altered paths (step 1706). If so, the 
routine resamples the altered paths (step 1708) and then constructs for each tooth an 

20 alternative path among the samples (step 1710). The routine continues in this manner 
until no collisions occur (step 1712). 

In some embodiments, as alluded to above, the software automatically 
computes the treatment path, based upon the IDDS and the FDDS. This is accomplished 
using a path scheduling algorithm which determines the rate at which each component, 

25 i.e., each tooth, moves along the path from the initial position to the final position. The 
path scheduling algorithm determines the treatment path while avoiding "round-tripping," 
i.e., while avoiding moving a tooth along a distance greater than absolutely necessary to 
straighten the teeth. Such motion is highly undesirable, and has potential negative effects 
on the patient. 

30 One implementation of the path scheduling algorithm attempts first to 

schedule or stage the movements of the teeth by constraining each tooth to the most linear 
treatment path between the initial and final positions. The algorithm then resorts to less 
direct routes to the final positions only if collisions will occur between teeth along the 
linear paths or if mandatory constraints will be violated. The algorithm applies one of the 
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path-generation processes described above, if necessary, to construct a path for which the 
intermediate treatment steps do not lie along a linear transformation curve between the 
initial and final positions. Alternatively, the algorithm schedules treatment paths by 
drawing upon a database of preferred treatments for exemplary tooth arrangements. This 
5 database can be constructed over time by observing various courses of treatment and 
identifying the treatment plans that prove most successful with each general class of 
initial tooth arrangements. The path scheduling algorithm can create several alternative 
paths and present each path graphically to the user. The algorithm provides as output the 
path selected by the user. 

10 In other implementations, the path scheduling algorithm utilizes a 

stochastic search technique to find an unobstructed path through a configuration space 
which describes possible treatment plans. One approach to scheduling motion between 
two user defined global key frames is described below. Scheduling over a time interval 
which includes intermediate key frames is accomplished by dividing the time interval into 

15 subintervals which do not include intermediate key frames, scheduling each of these 
intervals independently, and then concatenating the resulting schedules. 

Flow chart 120 in FIG. 8 A depicts a simplified path scheduling algorithm. 
As shown in FIG. 8 A, first step 122 involves construction of the "configuration space" 
description. A "configuration," in this context, refers to a given set of positions of all the 

20 teeth being considered for movement. Each of these positions may be described in 

multiple ways. In a common implementation, the positions are described by one affine 
transformation to specify change in location and one rotational transformation to specify 
the change in orientation of a tooth from its initial position to its final position. The 
intermediate positions of each tooth are described by a pair of numbers which specify 

25 how far to interpolate the location and orientation between the two endpoints. A 
"configuration" thus consists of two numbers for each tooth being moved, and the 
"configuration space" refers to the space of all such number pairs. Thus, the 
configuration space is a Cartesian space, any location in which can be interpreted as 
specifying the positions of all teeth. 

30 The affine transformation describing the movement of each tooth from its 

starting position to its ending position is decomposed into translational and rotational 
components; these transformations are independently interpolated with scalar parameters 
which are considered two dimensions of the configuration space. The entire 
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configuration space thus consists of two dimensions per moved tooth, all of which are 
treated equivalently during the subsequent search. 

The configuration space is made of "free space" and "obstructed space." 
"Free" configurations are those which represent valid, physically realizable positions of 
5 teeth, while "obstructed" configurations are those which do not. To determine whether a 
configuration is free or obstructed, a model is created for the positions of the teeth which 
the configuration describes. A collision detection algorithm is then applied to determine 
if any of the geometries describing the tooth surfaces intersect. If there are no 
obstructions, the space is considered free; otherwise it is obstructed. Suitable collision 

1 0 detection algorithms are discussed in more detail below. 

At step 124, a "visibility" function V(sj, S2) is defined which takes two 
vectors in the configuration space, "si" and M S2 f \ as input and returns a true or false 
boolean value. The visibility function returns a true value if and only if a straight line 
path connecting si and S2 passes entirely through a free and unobstructed region of the 

1 5 configuration space. One process for carrying out the visibility function is set forth in 

FIG. 8B. The visibility function is approximately computed by testing the teeth model for 
interferences at discretely sampled points along the line S1-S2. Techniques such as early 
termination on failure and choosing the order of sample points by recursively subdividing 
the interval to be tested, may be used to increase the efficiency of the visibility function. 

20 At step 126 of FIG. 8 A, a "children" function C(s) is defined whose input 

parameter, "s", is a vector in the configuration space and which returns a set of vectors 
"s c " in the configuration space. FIG. 8C depicts a simplified flow chart illustrating the 
steps followed for computing children function C(s). Each vector within set s c satisfies 
the property that V(s 5 s c ) is true and that each of its components are greater than or equal 

25 to the corresponding component of "s." This implies that any state represented by such a 
vector is reachable from "s" without encountering any interferences and without 
performing any motion which is not in the direction prescribed by treatment. Each vector 

of set "s c " is created by perturbing each component of "s" by some random, positive 

amount. The visibility function V(s s s c ) is then computed and "s" added to the set "s c " if 
30 the visibility function returns a true boolean value. Additionally, for each such vector 
generated, a pointer to its parent "s" is recorded for later use. 

After the configuration space has been defined, at step 128, path 
scheduling is performed between an initial state "Si n i t " and a final state "Sf ma r. FIG. 8D 
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depicts a flow chart for performing step 128 depicted in FIG. 8 A. As illustrated in FIG. 
8D, at step 128a 5 a set of states "W" is defined to initially contain only the initial state Si n i t . 
Next, at step 128b 3 the visibility function is invoked to determine if V(s, Sf ina i) is true for 
at least one state Si in W. If the visibility function returns a false boolean value, at step 
5 128c, the set of states "W" is replaced with the union of C(Si) for all Si in W. Steps 128b 
and 128c are repeated until V(si, Sf ma i) returns a true boolean value for any Si belonging 
to W. 

At step 128d ? for each Sj for which V(si, Sf ma i) is true, an unobstructed path 
Pi is constructed from Si to Si n u by following the parent pointers back to Sum. At step 1 28e, 

1 0 the path from Sjnj t to Sf ma i is then constructed by concatenating the paths Pi with the final 
step from Sj to s^ai- If there are multiple paths from Sinu to Sf ma i, the total length of each 
path is computed at step 128f. Finally, at step 128g, the path with the shortest length is 
then chosen as the final path. The length of the chosen path corresponds to the total time 
and stages required for a treatment plan. 

15 The resulting final path consists of a series of vectors, each of which 

represents a group of values of the interpolation parameters of the translational and 
rotational components of the transformations of the moving teeth. Taken together, these 
constitute a schedule of tooth movement which avoids tooth-to-tooth interferences. 

A collision or interference detection algorithm employed in one 

20 embodiment is based on the algorithm described in SIGGRAPH article, Stefan Gottschalk 
et al. (1996): "OBBTree: A Hierarchical Stincture for Rapid Interference Detection. " 
The contents of the SIGGRAPH article are herein incorporated by reference. 

The algorithm is centered around a recursive subdivision of the space 
occupied by an object, which is organized in a binary-tree like fashion. Triangles are 

25 used to represent the teeth in the DDS. Each node of the tree is referred to as an oriented 
bounding box (OBB) and contains a subset of triangles appearing in the node's parent. 
The children of a parent node contain between them all of the triangle data stored in the 
parent node. 

The bounding box of a node is oriented so it tightly fits around all of the 
30 triangles in that node. Leaf nodes in the tree ideally contain a single triangle, but can 
possibly contain more than one triangle. Detecting collisions between two objects 
involves determining if the OBB trees of the objects intersect. FIG. 9A sets forth a flow 
chart depicting a simplified version of a recursive collision test to check if a node "Nl" 
from a first object intersects with node "N2" of a second object. If the OBBs of the root 
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nodes of the trees overlap, the root's children are checked for overlap. The algorithm 
proceeds in a recursive fashion until the leaf nodes are reached. At this point, a robust 
triangle intersection routine is used to determine if the triangles at the leaves are involved 
in a collision. 

5 The collision detection technique described here provides several 

enhancements to the collision detection algorithm described in the SIGGRAPH article. 
For example, OBB trees can be built in a lazy fashion to save memory and time. This 
approach stems from the observation that some parts of the model will never be involved 
in a collision, and consequently the OBB tree for such parts of the model need not be 

1 0 computed. The OBB trees are expanded by splitting the internal nodes of the tree as 

necessary during the recursive collision determination algorithm, as depicted in FIG. 9B. 

Moreover, the triangles in the model which are not required for collision 
data may also be specifically excluded from consideration when building an OBB tree. 
As depicted in FIG. 9C, additional information is provided to the collision algorithm to 

15 specify objects in motion. Motion may be viewed at two levels. Objects may be 

conceptualized as "moving" in a global sense, or they may be conceptualized as "moving" 
relative to other objects. The additional information improves the time taken for the 
collision detection by avoiding recomputation of collision information between objects 
which are at rest relative to each other since the state of the collision between such objects 

20 does not change. 

FIG. 18 illustrates an alternative collision detection scheme, one which 
calculates a "collision buffer" oriented along a z-axis 1802 along which two teeth 1 804, 
1 806 lie. The collision buffer is calculated for each treatment step or at each position 
along a treatment path for which collision detection is required. To create the buffer, an 

25 x,y plane 1 808 is defined between the teeth 1 804, 1 806. The plane must be "neutral" with 
respect to the two teeth. Ideally, the neutral plane is positioned so that it does not 
intersect either tooth. If intersection with one or both teeth is inevitable, the neutral plane 
is oriented such that the teeth lie, as much as possible, on opposite sides of the plane. In 
other words, the neutral plane minimizes the amount of each tooth's surface area that lies 

30 on the same side of the plane as the other tooth. 

In the plane 1808 is a grid of discrete points, the resolution of which 
depends upon the required resolution for the collision detection routine. A typical high- 
resolution collision buffer includes a 400 x 400 grid; a typical low-resolution buffer 
includes a 20 x 20 grid. The z-axis 1 802 is defined by a line normal to the plane 1 808. 




The relative positions of the teeth 1 804, 1 806 are determined by 
calculating, for each of the points in the grid, the linear distance parallel to the z-axis 
1 802 between the plane 1 808 and the nearest surface of each tooth 1 804, 1 806. For 
example, at any given grid point (M,N), the plane 1808 and the nearest surface of the rear 
5 tooth 1 804 are separated by a distance represented by the value Zi(m,n)> while the plane 
1 808 and the nearest surface of the front tooth 1 806 are separated by a distance 
represented by the value Z 2 (m,n> If the collision buffer is defined such that the plane 1 808 
lies at z=0 and positive values of z lie toward the back tooth 1 804, then the teeth 1 804, 
1 806 collide when Zi (M ,n) ^ Z 2 (m,n) at any grid point (M,N) on the plane 1 808. 
10 FIG. 19 is a flow chart of a collision detection routine implementing this 

collision buffer scheme. The routine first receives data from one of the digital data sets 
indicating the positions of the surfaces of the teeth to be tested (step 1900). The routine 
then defines the neutral x,y-plane (step 1902) and creates the z-axis normal to the plane 
(step 1904). 

1 5 The routine then determines for the x,y-position of the first grid point on 

the plane the linear distance in the z-direction between the plane and the nearest surface 
of each tooth (step 1906). To detect a collision at that x,y-position, the routine determines 
whether the z-position of the nearest surface of the rear tooth is less than or equal to the z- 
position of the nearest surface of the front tooth (step 1908). If so, the routine creates an 

20 error message, for display to a user or for feedback to the path-generating program, 

indicating that a collision will occur (step 1910). The routine then determines whether it 
has tested all x,y-positions associated with grid points on the plane (step 1912) and, if not, 
repeats the steps above for each remaining grid point. The collision detection routine is 
performed for each pair of adjacent teeth in the patient's mouth at each treatment step. 

25 Incorporating a model of an orthodontic appliance . 

Above-mentioned U.S. application (attorney docket 

number 09943/004001) describes an appliance modeling system that implements 
techniques for modeling the interaction of the patient's teeth with orthodontic appliances 
designed to carry out the patient's treatment plan. Finite element analysis is used to 

30 determine the appliance configurations required to move the patient's teeth to the desired 
final positions along the prescribed treatment paths. In some situations, the appliance 
modeling system may determine that the desired tooth movement cannot be performed 
within constraints that are orthodontically acceptable or with an appliance that is 
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manufacturable. The appliance modeling system therefore may determine that a tooth 
attachment should be added to the model or that the treatment plan should be modified. 
In these situations, feedback from the appliance modeling system is used to modify the 
geometric tooth models and the treatment plan accordingly. 

5 DISPLAYING THE TREATMENT PLAN GRAPHICALLY 

The system may also incorporate and the user may at any point use a 
"movie" feature to show an animation of the movement from initial to target states. This 
is helpful for visualizing overall component movement throughout the treatment process. 
As described above, one suitable user interface for component 

1 0 identification is a three dimensional interactive graphical user interface (GUI). A three- 
dimensional GUI is also advantageous for component manipulation. Such an interface 
provides the treating professional or user with instant and visual interaction with the 
digital model components. The three-dimensional GUI provides advantages over 
interfaces that permit only simple low-level commands for directing the computer to 

15 manipulate a particular segment. In other words, a GUI adapted for manipulation is better 
in many ways than an interface that accepts directives, for example, only of the sort: 
"translate this component by 0.1 mm to the right." Such low-level commands are useful 
for fine-tuning, but, if they were the sole interface, the processes of component 
manipulation would become a tiresome and time-consuming interaction. 

20 Before or during the manipulation process, one or more tooth components 

may be augmented with template models of tooth roots. Manipulation of a tooth model 
augmented with a root template is useful, for example, in situations where impacting of 
teeth below the gumline is a concern. These template models could, for example, 
comprise a digitized representation of the patient's teeth x-rays. 

25 The software also allows for adding annotations to the data sets which can 

comprise text and/or the sequence number of the apparatus. The annotation is added as 
recessed text (i.e., it is 3-D geometry), so that it will appear on the printed positive model. 
If the annotation can be placed on a part of the mouth that will be covered by a 
repositioning appliance, but is unimportant for the tooth motion, the annotation may 

30 appear on the delivered repositioning appliance(s). 

The above-described component identification and component 
manipulation software is designed to operate at a sophistication commensurate with the 
operator's training level. For example, the component manipulation software can assist a 
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computer operator, lacking orthodontic training, by providing feedback regarding 
permissible and forbidden manipulations of the teeth. On the other hand, an orthodontist, 
having greater skill in intraoral physiology and teeth-moving dynamics, can simply use 
the component identification and manipulation software as a tool and disable or otherwise 
5 ignore the advice. 



with a client viewer application through which a treating clinician is able to view a 
patient's treatment plan and alter or comment on the plan. The client viewer application 
is implemented in a computer program installed locally on a client computer at the 

1 0 clinician's site. The viewer program downloads a data file from a remote host, such as a 
file transfer protocol (FTP) server maintained by the treatment plan designer, which can 
be accessed either through direct connection or through a computer network, such as the 
World Wide Web. The viewer program uses the downloaded file to present the treatment 
plan graphically to the clinician. The viewer program also can be used by the treatment 

1 5 plan designer at the host site to view images of a patient's teeth. 



treatment positions, including the IDDS and the FDDS, that define the treatment plan for 
the patient's teeth. The viewer program renders the IDDS or the FDDS to display an 
image of the patient's teeth at the initial and final positions. The viewer program can 
20 display an image of the teeth at their initial positions (initial image 2002) and the final 
tooth positions (final image 2004) concurrently. 



software in the remote host employs a "level-of-detail" technique to organize the 
download into data groups with progressively increasing levels of detail, as described 
25 below. The viewer program uses knowledge of orthodontic relevance to render less 

important areas of the image at a lower quality than it renders the more important areas. 
Use of these techniques reduces the time required to generate a single rendered image of 
the tooth models and the time required to display a rendered image on the screen after the 
download has begun. 

30 FIGS. 21 A and 21B illustrate the use of the "level-of-detail" technique by 

the download software in the remote host. The software transfers the data in several 
groups, each of which adds detail incrementally for the rendered image of the teeth. The 
first group typically includes just enough data to render a rough polygon representation of 
the patient's teeth. For example, if a tooth is treated as a cube having six faces, the tooth 



FIG. 20 is a screen shot of the graphical user interface 2000 associated 



The data downloaded by the viewer program contains a fixed subset of key 



Because the data file contains a large amount of data, the download 
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can be rendered quickly as a diamond 2100 having six points 2102a-f, one lying in each 
face of the cube (FIG. 21 A). The download software begins the download by delivering a 
few points for each tooth, which the interface program uses to render polygon 
representations of the teeth immediately. 
5 The download software then delivers a second data group that adds 

additional detail to the rendered images of the teeth. This group typically adds points that 
allow a spheroid representation 2106 of the teeth (FIG. 2 IB). As the download continues, 
the software delivers additional groups of data, each adding a level of detail to the 
rendered image of the teeth, until the teeth are fully rendered. 

1 0 The download software also improves download and rendering speed by 

identifying and withholding data that is not critical to forming a rendered image of the 
teeth. This includes data for tooth surfaces obscured by other teeth or by tissue. The 
software applies rules based on common orthodontic structure to determine which data is 
downloaded and which is withheld. Withholding data in this manner reduces the size of 

1 5 the downloaded file and therefore reduces the number of data points that the interface 
program must take into account when rendering the initial and final images. 

The viewer program also improves rendering speed by reducing the 
amount of data rendered. Like the download software, the viewer program applies rules 
of orthodontic relevance to determine which areas of the image can be rendered at lower 

20 quality. For example, the treating clinician usually does not want to view gum tissue in 
detail, so the viewer program renders the gums at low resolution as smooth surfaces, 
ignoring data that preserves the texture of the gums. Typically, the viewer program 
renders the less important areas at lower resolution before rendering the more important 
areas at higher resolution. The clinician can request high resolution rendering of the 

25 entire image. 

As shown in FIG. 20 and discussed above, the viewer program displays an 
initial image 2002 of the teeth and, if requested by the clinician, a final image 2004 of the 
teeth as they will appear after treatment. The clinician can rotate the images in three 
dimensions to view the various tooth surfaces, and the clinician can snap the image to any 
30 of several predefined viewing angles. These viewing angles include the standard front, 
back, top, bottom and side views, as well as orthodontic-specific viewing angles, such as 
the lingual, buccal, facial, occlusal, and incisal views. 

The viewer program also includes an animation routine that provides a 
series of images showing the positions of the teeth at each intermediate step along the 
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treatment path. The clinician controls the animation routine through a VCR metaphor, 
which provides control buttons similar to those on a conventional video cassette recorder. 
In particular, the VCR metaphor includes a "play" button 2006 that, when selected, 
causes the animation routine to step through all of the images along the treatment path. A 
5 slide bar 2008 moves horizontally a predetermined distance with each successive image 
displayed. Each position of the slide bar 2008 and each image in the series corresponds 
to one of the intermediate treatment steps described above. 

The VCR metaphor also includes a "step forward" button 2010 and a "step 
back" button 2012, which allow the clinician to step forward or backward through the 

1 0 series of images, one key frame or treatment step at a time, as well as a "fast forward" 
button 2014 and a "fast back" button 201 6, which allow the clinician to jump 
immediately to the final image 2004 or initial image 2002, respectively. The clinician 
also can step immediately to any image in the series by positioning the slide bar 2008 at 
the appropriate location. 

15 As described above, the viewer program receives a fixed subset of key 

positions, including the IDDS and the FDDS, from the remote host. From this data, the 
animation routine derives the transformation curves required to display the teeth at the 
intermediate treatment steps, using any of a variety of mathematical techniques. One 
technique is by invoking the path-generation program described above. In this situation, 

20 the viewer program includes the path-generation program code. The animation routine 

invokes this code either when the downloaded key positions are first received or when the 
user invokes the animation routine. 

The viewer program allows the clinician to alter the rendered image by 
manipulating the image graphically. For example, the clinician can reposition an 

25 individual tooth by using a mouse to click and drag or rotate the tooth to a desired 
position. In some implementations, repositioning an individual tooth alters only the 
rendered image; in other implementations, repositioning a tooth in this manner modifies 
the underlying data set. In the latter situation, the viewer program performs collision 
detection to determine whether the attempted alteration is valid and, if not, notifies the 

30 clinician immediately. Alternatively, the viewer program modifies the underlying data 
set and then uploads the altered data set to the remote host, which performs the collision 
detection algorithm. The clinician also can provide textual feedback to the remote host 
through a dialog box 2018 in the interface display 2000. Text entered into the dialog box 
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2018 is stored as a text object and later uploaded to the remote host or, alternatively, is 
delivered to the remote host immediately via an existing connection. 

The viewer program optionally allows the clinician to isolate the image of 
a particular tooth and view the tooth apart from the other teeth. The clinician also can 
5 change the color of an individual tooth or group of teeth in a single rendered image or 

across the series of images. These features give the clinician a better understanding of the 
behavior of individual teeth during the course of treatment. 

Another feature of the viewer program allows the clinician to receive 
information about a specific tooth or a specific part of the model upon command, e.g., by 
1 0 selecting the area of interest with a mouse. The types of information available include 
tooth type, distance between adjacent teeth, and forces (magnitudes and directions) 
exerted on the teeth by the aligner or by other teeth. Finite element analysis techniques 
are used to calculate the forces exerted on the teeth. The clinician also can request 
graphical displays of certain information, such as a plot of the forces exerted on a tooth 
1 5 throughout the course of treatment or a chart showing the movements that a tooth will 
make between steps on the treatment path. The viewer program also optionally includes 
"virtual calipers," a graphical tool that allows the clinician to select two points on the 
rendered image and receive a display indicating the distance between the points. 

FABRICATING THE ALIGNERS 
20 Once the intermediate and final data sets have been created, the appliances 

may be fabricated as illustrated in FIG. 10. Common fabrication methods employ a rapid 
prototyping device 200 such as a stereolithography machine. A particularly suitable rapid 
prototyping machine is Model SLA-250/50 available from 3D System, Valencia, 
California. The rapid prototyping machine 200 selectively hardens a liquid or other non- 
25 hardened resin into a three-dimensional structure which can be separated from the 
remaining non-hardened resin, washed, and used either directly as the appliance or 
indirectly as a mold for producing the appliance. The prototyping machine 200 receives 
the individual digital data sets and produces one structure corresponding to each of the 
desired appliances. Generally, because the rapid prototyping machine 200 may utilize a 
30 resin having non-optimum mechanical properties and which may not be generally 

acceptable for patient use, the prototyping machine typically is used to produce molds 
which are, in effect, positive tooth models of each successive stage of the treatment. 
After the positive models are prepared, a conventional pressure or vacuum molding 
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machine is used to produce the appliances from a more suitable material, such as 0.03 
inch thermal forming dental material, available from Tru-Tain Plastics, Rochester, 
Minnesota 55902. Suitable pressure molding equipment is available under the trade name 
BIOSTAR from Great Lakes Orthodontics, Ltd., Tonawanda, New York 14150. The 
5 molding machine 250 produces each of the appliances directly from the positive tooth 
model and the desired material. Suitable vacuum molding machines are available from 
Raintree Essix, Inc. 



professional all at one time. The appliances are marked in some manner, typically by 
10 sequential numbering directly on the appliances or on tags, pouches, or other items which 
are affixed to or which enclose each appliance, to indicate their order of use. Optionally, 
written instructions may accompany the system which set forth that the patient is to wear 
the individual appliances in the order marked on the appliances or elsewhere in the 
packaging. Use of the appliances in such a manner will reposition the patient's teeth 
1 5 progressively toward the final tooth arrangement. 



the treating clinician may wish to evaluate the patient's progress during the course of 
treatment. The system can also do this automatically, starting from the newly-measured 
in-course dentition. If the patient's teeth do not progress as planned, the clinician can 

20 revise the treatment plan as necessary to bring the patient's treatment back on course or to 
design an alternative treatment plan. The clinician may provide comments, oral or 
written, for use in revising the treatment plan. The clinician also can form another set of 
plaster castings of the patient's teeth for digital imaging and manipulation. The clinician 
may wish to limit initial aligner production to only a few aligners, delaying production on 

25 subsequent aligners until the patient's progress has been evaluated. 



may be used to develop orthodontic treatment plans. The data processing system 300 
typically includes at least one processor 302 which communicates with a number of 
peripheral devices via bus subsystem 304. These peripheral devices typically include a 
30 storage subsystem 306 (memory subsystem 308 and file storage subsystem 314), a set of 
user interface input and output devices 318, and an interface to outside networks 3 1 6, 
including the public switched telephone network. This interface is shown schematically 
as "Modems and Network Interface" block 316, and is coupled to corresponding interface 
devices in other data processing systems via communication network interface 324. Data 



After production, the appliances can be supplied to the treating 



Because a patient's teeth may respond differently than originally expected, 



FIG. 1 1 is a simplified block diagram of a data processing system 300 that 
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processing system 300 could be a terminal or a low-end personal computer or a high-end 
personal computer, workstation or mainframe. 

The user interface input devices typically include a keyboard and may 
further include a pointing device and a scanner. The pointing device may be an indirect 
5 pointing device such as a mouse, trackball, touchpad, or graphics tablet, or a direct 
pointing device such as a touchscreen incorporated into the display, or a three 
dimensional pointing device, such as the gyroscopic pointing device described in U.S. 
Patent 5,440,326, other types of user interface input devices, such as voice recognition 
systems, can also be used. 

10 User interface output devices typically include a printer and a display 

subsystem, which includes a display controller and a display device coupled to the 
controller. The display device may be a cathode ray tube (CRT), a flat-panel device such 
as a liquid crystal display (LCD), or a projection device. The display subsystem may also 
provide non-visual display such as audio output. 

1 5 Storage subsystem 306 maintains the basic required programming and data 

constructs. The program modules discussed above are typically stored in storage 
subsystem 306. Storage subsystem 306 typically comprises memory subsystem 308 and 
file storage subsystem 314. 

Memory subsystem 308 typically includes a number of memories 

20 including a main random access memory (RAM) 3 1 0 for storage of instructions and data 
during program execution and a read only memory (ROM) 3 12 in which fixed 
instructions are stored. In the case of Macintosh-compatible personal computers the 
ROM would include portions of the operating system; in the case of IBM-compatible 
personal computers, this would include the BIOS (basic input/output system). 

25 File storage subsystem 314 provides persistent (non- volatile) storage for 

program and data files, and typically includes at least one hard disk drive and at least one 
floppy disk drive (with associated removable media). There may also be other devices 
such as a CD-ROM drive and optical drives (all with their associated removable media). 
Additionally, the system may include drives of the type with removable media cartridges. 

30 The removable media cartridges may, for example be hard disk cartridges, such as those 
marketed by Syquest and others, and flexible disk cartridges, such as those marketed by 
Iomega. One or more of the drives may be located at a remote location, such as in a 
server on a local area network or at a site on the Internets World Wide Web. 
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In this context, the term "bus subsystem" is used genetically so as to 



include any mechanism for letting the various components and subsystems communicate 
with each other as intended. With the exception of the input devices and the display, the 
other components need not be at the same physical location. Thus, for example, portions 
5 of the file storage system could be connected via various local-area or wide-area network 
media, including telephone lines. Similarly, the input devices and display need not be at 
the same location as the processor, although it is anticipated that personal computers and 
workstations typically will be used. 



1 0 system has a number of buses such as a local bus and one or more expansion buses (e.g., 
ADB, SCSI, ISA, EISA, MCA, NuBus, or PCI), as well as serial and parallel ports. 
Network connections are usually established through a device such as a network adapter 
on one of these expansion buses or a modem on a serial port. The client computer may be 
a desktop system or a portable system. 

15 Scanner 320 is responsible for scanning casts of the patient ! s teeth obtained 

either from the patient or from an orthodontist and providing the scanned digital data set 
information to data processing system 300 for further processing. In a distributed 
environment, scanner 320 may be located at a remote location and communicate scanned 
digital data set information to data processing system 300 via network interface 324. 

20 Fabrication machine 322 fabricates dental appliances based on 

intermediate and final data set information received from data processing system 300. In 
a distributed environment, fabrication machine 322 may be located at a remote location 
and receive data set information from data processing system 300 via network interface 



Other embodiments are within the scope of the following claims. For example, the three- 
dimensional scanning techniques described above may be used to analyze material 
characteristics, such as shrinkage and expansion, of the materials that form the tooth 
castings and the aligners. Also, the 3D tooth models and the graphical interface described 
30 above may be used to assist clinicians that treat patients with conventional braces or other 
conventional orthodontic appliances, in which case the constraints applied to tooth 
movement would be modified accordingly. Moreover, the tooth models may be posted 
on a hypertext transfer protocol (http) web site for limited access by the corresponding 
patients and treating clinicians. 



Bus subsystem 304 is shown schematically as a single bus, but a typical 



324. 



25 



The invention has been described in terms of particular embodiments. 
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